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Abstract 

We consider the effect of differing coefficients of static and dynamic friction 
coefficients on the behaviour of contacts involving microslip. The classic solu¬ 
tions of Cattaneo and Mindlin are unchanged if the transition in coefficients is 
abrupt, but if it occurs over some small slip distance, the solution has some 
mathematical similarities with those governing the normal tractions in adhe¬ 
sive contact problems. In particular, if the transition to dynamic slip occurs 
over a sufficiently small area, we can identify a ‘ JKR’ approximation, where the 
transition region is condensed to a line. A local singularity in shear traction 
is then predicted, with a stress-intensity factor that is proportional to the the 
square root of the local contact pressure and to a certain integral of the friction 
coefficient-slip distance relation. We can also define an equivalent of the ‘small- 
scale yielding’ criterion, which enables us to assess when the singular solution 
provides a good approximation. One consequence of the results is that the 
static coefficient of friction determined from force measurements in experiments 
is significantly smaller than the value that holds at the microscale. 

1 Introduction 

If a deformable structure with frictional interfaces is subjected to loads that 
are insufficient to cause gross slip (sliding), the deformation of the components 
generally permits some local regions of ‘microslip’ at the nominally stuck contact 
interfaces. When the loading is periodic, these regions contribute to the energy 
dissipation in the structure and hence influence the dynamic behaviour mm. 
Also, cyclic microslip can eventually lead to the initiation and propagation of 
fretting fatigue cracks [5]. 

Most of the extensive literature on problems involving microslip assumes 
that Coulomb’s friction law applies — i.e. 

u 

q = -/pttt; m/o (i) 

M 

l<?l < fp ; u = o, (2) 

where q is the frictional (tangential) traction, p is the contact pressure, u is 
the local microslip velocity, and / is the coefficient of friction. In particular, it 
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is usually assumed that the same coefficient / governs both the slip and stick 
regions. 

By contrast, dynamicists and tribologists often make a distinction between 
static and dynamic friction [2], so that equations m are replaced by 

u 

Q = —fdP 7tt ; ft 7^0 (3) 

M 

\q\ < fsp; ft = 0, (4) 

where f s , fd are the static and dynamic friction coefficients respectively. In par¬ 
ticular, if f s > fd, this friction law provides a mechanism for ‘stick-slip’ frictional 
vibrations [5]. Numerous experimental investigations have shown differences be¬ 
tween static and sliding friction (e.g. 0 ). These differences are generally small 
for dry metals [7j, but can be substantial for earthquake fault mechanics, where 
ratios as high as ten between the coefficients have been reported [8]. Rice 0 
characterizes such interfaces as ‘strong but brittle’. 

A higher coefficient of static friction can to some extent be explained by 
noting that the formation of adhesive bonds, which forms the basis of Bowden 
and Tabor’s friction theory m, will be enhanced by diffusion if asperities remain 
in contact for some period of time. Similar arguments can be used to justify the 
‘rate-state’ friction model mm 

In this paper, we shall examine the effect of introducing a higher coefficient 
of static friction on problems involving microslip. In the interests of simplicity, 
we shall restrict attention to cases where Dundurs’ parameter /3 = 0 ([13] p. 
110), so there is no coupling between normal and tangential loading, and the 
contact pressure can be determined without reference to the friction law. Also, 
we shall illustrate our ideas in the context of the two-dimensional Hertz problem, 
since this is susceptible to simple analytical solutions, but extension to other 
two-dimensional cases, and to the axisymmetric Hertz problem is routine. 


2 Evolution of frictional traction distributions 

Cattaneo [14] and later Mindlin m considered the case where two elastic 
bodies are first pressed together by a normal force P, which is then held 
constant whilst a monotonically increasing unidirectional force Q x is applied. 
The profile of the bodies was characterized by a quadratic initial gap function 
go{x,y ) = Ax 2 + By 2 , so that the normal loading phase is defined by the clas¬ 
sical Hertz theory. Cattaneo and Mindlin then showed that, subject to a small 
approximation associated with the local slip direction m, the shear traction 
distribution has the form 

Qx{x,y) = f [p(x, y ) -p*(x, y)] , (5) 

where p(x,y) is the contact pressure and p*(x,y ) is the contact pressure that 
would be developed at some smaller normal force P* given by 

P* =P-9a . (6) 
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Ciavarella m and Jager [IS] have since shown that this form of superposition 
is exact for any initial gap function go(x) in the two dimensional case, and that 
it is a good approximation in the general three-dimensional case nn. 


3 Static and dynamic friction 

Now consider the case where f s > f d and the loading scenario is the same as in 
the Cattaneo-Mindlin problem. We assume the existence of a slip zone in which 
q x {x,y) = fdp(x,y), so we write the complete shear traction distribution as 

Qx(x,y) = f d p{x,y) - q* x (x,y) , (7) 

where q*(x,y) is a corrective distribution to be determined from the condition 
that the slip displacement (i.e. the relative tangential displacement) is zero in 
the stick area *4 s ti c k- Conditions (13171) require that q*(x,y) be non-zero only 
in -Astickj and hence the stick condition defines a well-posed boundary-value 
problem for q*(x,y). The inequality condition © precludes singularities in 
the shear tractions, and this imposes uniqueness on the solution for any given 
Qx < fd P■ It is clear that the original Cattaneo-Mindlin solution © with 
f = fd satisfies these conditions, including the inequality, since in _4 s tick, this 
would give q x {x,y) < fdp(x,y) < f s p(x,y). 


4 Dependence on slip distance 


The discussion so far is predicated on the assumption that as soon as stick is 
‘broken’ there is an immediate transition to the dynamic coefficient fd, but in 
practice we might expect a more continuous transition as slip occurs. We shall 
therefore examine the consequences of a friction law in which the coefficient of 
friction is a continuous and monotonic function f(u) of the slip displacement u, 
such that 

/(0) = f s and /(it) -)• f d ; u oo . (8) 


Such a law can be regarded as a special case of the rate-state law muni and 
is also related to the the shear failure law proposed by Abercrombie and Rice 
[20] , Applications of similar laws to fault mechanics are discussed by Ben Zion 

Ell¬ 


in general, solutions of the corresponding contact problem will then require 
numerical solution, but it is instructive to consider some simple cases analyti¬ 
cally. In particular, we shall consider the two-dimensional case where the bodies 
comprise a cylinder of radius R and a half space, so the contact pressure is given 
by 


p(x) 


E * y/ a 2 — x 2 7 tE* a 2 

2 R ’ = 4R 


( 9 ) 


where a is the semi-width of the contact area —a < x < a, and E is the 
composite elastic modulus 113]. 
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We anticipate the existence of two symmetric slip regions —a<x<—c and 
c < x < a in which the slip displacement increases monotonically away from the 
stick-slip boundaries x = ±c. Two limiting cases can also be identified. If f(u) 
is a rather slowly decaying function of u, the friction coefficient will be close to 
f s throughout the slip regions and the solution will approximate the constant 
coefficient case with / = f s . At the other limit, if a very small amount of slip 
displacement is required to precipitate the change in coefficient, most of the slip 
area will be at or near fd , but we must still allow for the existence of ‘transition’ 
regions c < \x\ < b m which f > fd- 

The exact form of the function f(u) is not critical, but it is convenient to 
define a quantity W with the dimensions of surface energy through the relation 

p OO 

w= (f(u) - fd)pdu , (10) 

Jo 

which is equivalent to the shear fracture energy defined by Abercrombie and 
Rice [2D]. The contact pressure p will generally vary in the transition region, 
but if this is sufficiently short for p to be regarded as uniform, we can also 
define a length scale A characterizing the amount of slip needed to transition 
to dynamic friction, such that 


A = 


W 

(fs - fd)p 


1 


(fs - fd) 


(f(u) - fd) du 


( 11 ) 


Rabinowicz [B] conducted some simple but elegant experiments to determine 
f s , fd and A for metals, his result^ being presented in Table 1. 


Materials 

fs 

fd 

A (/im) 

copper/mild steel 

0.46 

0.31 

1 

lead/mild steel 

0.72 

0.47 

3 

mild steel/copper 

0.54 

0.39 

0.9 

mild steel/titanium 

0.63 

0.45 

6 

mild steel/zinc 

0.65 

0.47 

2 


Table 1: Friction coefficients and slip length A for some metal combinations, 
from [6] 

A special case satisfying equations mm is the step function f = f s — (f s — 
fd)H(u — A), where H(-) is the Heaviside step function. The perceptive reader 
will notice a similarity here to Maugis’ approximate formulation of the normal 
adhesive contact problem [22] . where the adhesion law is also represented by a 
step function and the outer boundary of the adhered region is determined from 
the condition that the separation there is equal to a critical value. Indeed we 
shall see that there are significant mathematical analogies between the present 
problem and adhesive problems. 

Ht is difficult to explain why different results might be obtained by simply interchanging 
the materials in the mild steel/copper case, but the difference is arguable within the range of 
likely experimental variance. 
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4.1 A double-Cattaneo-Mindlin solution 


The present problem could be formulated using a step function for f(u), but a 
simpler mathematical approximation can be obtained by adapting the ‘double- 
Hertz’ concept of Greenwood and Johnson [ 22 ] • We first note that the Cattaneo- 
Mindlin traction distribution q x (x) = q{x,a,c),Q x = Q(a,c), where 

/ 2 2 \ 

q(x , a, c) = \J a 2 — x 2 — \/c 2 — x 2 ; Q{a,c) =— — (12) 


produces slip displacements u x (x), such that 


du x _ , x n 

— = v{x,a,c) = 0; 


2\/x 2 — c 2 

F 


— C < X < c 


; c < |x| < a 


(13) 

(14) 


m p. 214), where the square roots in (1121) are to be interpreted as zero in any 
region where their respective arguments are negative. 

We next approximate the solution to the frictional problem as 


q x {x) 


E* fd q(x, a, c) 
2 R 


Cq{x , b , c) ; Q x 


E*fdQ(a,c) 
2 R 


CQ(b,c) , (15) 


where c < b < a. The corresponding slip displacements will then satisfy 



E* f d v{x, a, c) 
2 R 


Cv(x, 6, c) , 


(16) 


and this is zero in — c < x < c from m, showing that the stick condition can 
be satisfied by an appropriate rigid-body translation. 

The shear tractions in b < |x| < a are 


qx(x) 


E* f d V a, 2 - x 2 
2 R 


fdp{x) , 


(17) 


and hence satisfy the slip condition at f = fd, since the other square-root terms 
make no contribution in this range. In c < |x| < b, the shear tractions are 


qx{x) = fdP{x) + C\/b 2 - x 2 , 


(18) 


and we can choose the constant C so as to ensure that q x (c) = f s p(c ), giving 


E* (f s — fd) 

2 R \b 2 - c 2 


(19) 


and 


q x {x) 


E* fdq(x, a, c) 
2 R 


E* {fs — fd) 

2 R 



q{x,b,c) . 


( 20 ) 
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With this choice, the effective local coefficient of friction f = q x /p will decrease 
nronotonically from f s to fd in c < \x\ < b. 

The final step is to determine the unknown radii c, b from the equilibrium 
condition (1151) 9. and from m which we can write as 

/ b j 

[qx{x) - fdp(x)} dx . (21) 


In c < |x| < b , we have 

(/-</«< 22 > 

from (|14I19I) . Using this expression and (jTSll in (1211) and evaluating the integral, 
we obtain 


W = - 


E b(f s - f d ) a? -c 


6 R 2 


b 2 — c 2 


fd + (fs ~ fd) 


b 2 — c 2 


x [(b 2 + c 2 )E(k) - 2 c 2 K(k)\ 


where 


k 2 = 1 - 


b 2 


and 


n 2 

K(k) = / 

Jo 


d.0 


o \J\ — k 2 sin 2 9 


; E(k) = [ ' Vl-k 2 
Jo 


sin 2 0 dd 


(23) 

(24) 

(25) 


are the complete elliptic integrals of the first and second kind respectively. The 
equilibrium condition is obtained from (1121151191) as 


Qx — 


7 lE 

~AR. L 


fd(a 2 - c 2 ) + (f s - fd)\J(b 2 — c 2 )(a 2 - c 2 ) 


(26) 


If Q X ,W are given, 
c, b. 


provide two equations for the two unknown radii 


4.2 The ‘JKR’ limit 


If the transition from f s to fd occurs over a sufficiently small region, we can 
obtain a limiting solution analogous to the JKR solution of normal adhesion 
problems. We write b = c + <5, where iCc, in which case (1231) can be approxi¬ 
mated as 


W 


^E* (f s — f d ) 2 (a 2 — c 2 )5 
16 R 2 


implying 6 


16 R 2 W 


^E*(f s — f d ) 2 (a 2 — c 2 ) ' 
(27) 
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Also, the second term in q x {x) in equation (fl5l) can be approximated as 


Cq{x,b,c) « Cq(x,c,c) + C6—(x,c,c) = _ ^ . (28) 

Applying the same approximation to equations (11911261) and substituting for S 
from (I27[) . we obtain 


/ \ _ E fd q(x, a, c) 
qx{x) «-—-h 


2 R 


and 


Qx 

JdP ttE* f d a 2 


A R Q X _ c 2 

* j . _0 ~ 1 ,,2 


2WE c 

(29) 


7 — x z ) 


4 R \2Wc 

(30) 

fd a 2 V 7 tE* ' 


Equation (1291) defines a locally singular field, implying the existence of a 
mode II stress-intensity factor 


it'll = lim q x (x)y/2Tr(c — x ) = V2VFE* , 


(31) 


which is exactly analogous with the mode I stress intensity factor K\ = \J 2A 7 IS* 
in normal adhesion problems in the JKR limit, where A 7 is the interface energy. 

In an impressive series of experiments, Svetlizky and Fineberg [24] have 
observed frictional slip progressing by the relatively slow propagation of slip 
zones behind which the shear tractions approximate a square-root singularity. 
The strength of this singularity is approximately constant, indicating a well- 
defined value of fracture energy W, but they suggest it may depend on the 
local pressure, as a result of the area of actual contact being approximately 
proportional to pressure. 

Ciavarclla [23] presented solutions of contact problems with a mode II stress- 
intensity factor around the stick-slip boundary, motivated by Fineberg’s obser¬ 
vations. The present analysis shows that such an effect can be generated by a 
slip-dependent friction law of the form (0 and provides a rationale for deter¬ 
mining an appropriate value of A'n- In particular, we notice from (1311) that the 
stress-intensity factor depends only on the composite modulus and W, and is 
otherwise independent of the details of the contact problem. Since ex hypothesis 
the transition is assumed to occur over a small region (of width S) in the contact 
area, we can assume that the contact pressure p is uniform in this region, and 
hence use the form CD for W. This leads to a stress-intensity factor 

Ku = y/2E*(f a - f d ) P A , (32) 

which varies with yjp and is equivalent to the ‘pressure-dependent toughness’ 
criterion of [25]. 
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Using CP to recast equations CUES) in terms of A, we have 


, ^ E* f d q(x, a, c) * I (f s - f d )AcVa 2 - c 2 

'• (I) “ — 2 r — + E V „n (< ? - x 2 ) — 

Qx _ 1 C 2 4 / (/ s - fd)RAcVa 2 - c 2 

/d-P a 2 f d a 2 V 7r 


(33) 

(34) 


4.3 Small-scale transition zone 


Equation (13^1) implies that at a sufficiently small distance s from the stick bound¬ 
ary, the frictional tractions have the singular form 


Qx(s) « f d p + 


E*(f„-f d )pA 


(35) 


However, this expression violates the stick condition dj) in the region 0 < s < Sq, 
where 


< E*(f s -f d )pA 

TTS 0 


= (fs - fd)p or s 0 = 


E* A 


7r(/ S - fd)p ' 


(36) 


An analogous situation is encountered in elastic-plastic fracture mechanics, 
where the ‘small-scale yielding’ criterion is used to determine whether the fields 
far outside the yield zone can reasonably be described by the elastic solution 
P5] . In the present case, the singular solution can be expected to give good 
results everywhere except very close to x = c, provided so -C c. 

This criterion depends on c and hence on Q x , but a rough estimate of the 
applicability of the JKR solution in the present problem can be obtained by 
using p( 0), a for p, c respectively, defining the modified criterion 


A = 


RA 


ifs - fd)a 2 


< 1 


(37) 


4.4 More general two-dimensional problems 

We have analyzed the two-dimensional Hertzian problem in detail because the 
resulting expressions are algebraically straightforward, enabling the fundamen¬ 
tal structure of the solution to be exposed. However, the same method can be 
applied to any two-dimensional problem involving a single symmetric contact 
area. We simply replace equation (E|) by 

q(x, a, c) = p(x, a) — p(x, c) ; Q(a,c) = P(a) - P(c) , (38) 

where p(x , a) is the normal contact pressure when the contact area is defined 
by —a < x < a, and P(a) is the corresponding normal force. We know from 
Ciavarclla m and Jager [18] that this will satisfy equation so the traction 
distribution 

q x (x) = f d q(x,a,c) - Cq(x,b,c) (39) 
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will satisfy the stick conditions in —c < x < c and the dynamic slip conditions 
in b < \x\ < a. The rest of the solution can then be completed as in §4.1. 

If the length scale so hr (1551) is sufficiently small to justify the JKR approx¬ 
imation, the second term will take the universal form (1281) . so the solution can 
be written down as the superposition of a conventional Cattaneo-Mindlin solu¬ 
tion with coefficient of friction fd and equation (1281) . In this context, it may be 
helpful to note that the limiting expression for Q(b, c) is 

Q(b,c)= f Cq(x , 6, c)dx — > V2nWE*c , (40) 

J —c 

so the total tangential force is 

Q x = fd [P(a ) - P(c)] + V2nWE*c . (41) 

Since Q x will usually be prescribed, this provides an equation from which c can 
be determined as a function of Q x . 

5 Finite element results 

The double Cattaneo-Mindlin solution is approximate in the sense that we are 
able to match a specific value of the fracture energy W or (equivalently) the 
length scale A, but the exact form of the function f(u) cannot be prescribed. 
The implied form of this function depends on the dimensionless ratios 6/c, a/c, 
some representative curves being shown as dashed and dotted lines in Fig. 1. 



u 

A 


Fig. 1: The friction coefficient function f(u) implied by the double Cattaneo- 
Mindlin solution for a/c = 1.6, b/c = 1.5 (— ■ — ■ —), a/c = 8.0, b/c = 4.5 (■ ■ ■ ■), 
a/c = 8.0, b/c = 1.2 (-). 

To assess the effect of this approximation, we constructed a finite element 
solution of the problem, as an extension of the ’’verification manual” VM272 
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example in Ansys 15 m, which in turn is based on the method of Yang et 
al. [28] and an example given therein which compares satisfactorily with the 
analytical Cattaneo-Mindlin solution. It is based on a mortar formulation of 
the contact which is able to deal with nonconforming discretizations across 
boundaries and large sliding which is more than adequate for our problem. In 
[25] , several examples and comparisons are made to show that this method has 
an optimal convergence rate and robustness with respect to other approaches. 

The example considers two parallel linear elastic half cylinders of radius R and 
pressed by a small distributed pressure on the diameter. A tangential pressure is 
then applied to cause friction at the contact interface, while the top of the upper 
cylinder is constrained from rotating. The bottom of the lower cylinder is fixed 
in all directions. The standard input listing available in ANSYS is adequate for 
many problems, but two minor changes made in the present case were:- 

1. We used quadratic PLANE183-CONTA172 instead of linear elements PLANE182- 
CONTA171, and we modified the mesh parametrically keeping the same 

ratio of elements, in order to improve marginally the accuracy of the re¬ 
sults. For the figures reported in the paper we divided every element edge 
by 3 which brings the total number of elements to about 45000, but still 
permits a solution of an entire curve of loading in less than a minute. 

2. We did not use the ANSYS variant of the friction law with just static 
and dynamic coefficients, since this does not permit a dependence on slip 
displacement. Instead., we defined a table of friction coefficients in terms 
of slip displacement. 

Fig. 2 compares the shear traction q x (x) from equation (I2U1) for Q x = 

0.8 fdP, f s = 0.15, fd = 0.1, A = 0.05, with finite element results using the ramp 
(linear) function for f(u) from Fig. 1. The agreement is clearly extremely good. 

Also shown on this figure are the conventional Cattaneo-Mindlin prediction 
(equivalent to taking A = 0) and the JKR approximation (l33l) . The latter gives 
good predictions everywhere except in the transition region, where of course the 
predicted singular stress is unphysical. 
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Fig. 2: Finite element results (•••) for the shear traction distribution q x (x) for 
Q x = 0.8 fdP,fs = 0.15, fd = 0.1 and A = 0.05: ““ Double Cattaneo-Mindlin 

solution (EU1) . - Conventional Cattaneo-Mindlin solution (0 with / = fd, 

- ‘JKR’ approximation of equation (l33l) . 

Fig. 3 shows a similar comparison for a larger value of A, so that the 
transition extends over a larger radius. In this figure, we compare equation 
(1201) with finite element solutions using the ramp function and the step function 
respectively from Fig. 1. This figure shows that the traction distribution is 
relatively insensitive the the form of the function f{u ) for given values of (f s — fd) 
and A, and hence that equation (l20l) can be expected to give good results for 
most practical slip-weakening laws. 



Fig. 3: Effect of the function f(u) on the traction distribution q x {x): ••• ramp 

(finite element), - step (finite element), — equation (|20|) . Q x = 0.9 fdP, 

f. = 0.15, f d = 0.1, A = 0.277. 

6 Discussion 

The principal new result from this analysis is that fracture mechanics concepts 
are introduced into the microslip problem, even when the friction law is merely 
an extension of the Coulomb law allowing differing static and dynamic coeffi¬ 
cients. In particular, if the coefficient of friction varies with slip dispacement 
over a relatively short slip distance A, we can determine a critical stress in¬ 
tensity factor or fracture toughness (132)) that depends only on the static and 
dynamic coefficients, the form of the slip-weakening law, the composite modulus 
and the local contact pressure. 

Equation (l34l) defines the relation between tangential force Q x and the semi¬ 
length c of the stick area in the JKR limit, which is appropriate if the small-scale 
transition criterion Sq -C c is satisfied. It is plotted in Fig. 4 for several values 
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of the dimensionless parameter 


4- = 



(42) 



a 

Fig. 4 : The tangential force Q x as a function of the radius c of the stick zone 
(JKR limit). 

All the curves except the limiting case 4' = 0 exhibit a maximum Q x = Q max 
at a non-zero value of c, implying that under tangential force control, the system 
would jump unstably to full sliding once this maximum is reached. The unstable 
range is shown dotted in Fig. 4. 

Similar plots were made for the double Cattaneo-Mindlin solution, using 
equations (l26l l23l) with W = (/ s — f p ) A. Fig. 5 compares the resulting curves 
for ip = 0.1 and A = 0.025,0.4 with the JKR solution. Notice that changing 
A at constant ip implies a change in the friction coefficient ratio f s /fd■ The 
truncation in these curves near c = a occurs because the outer boundary b of 
the transition region cannot exceed the boundary a of the contact area. When 
b = a, the double Cattaneo-Mindlin solution reduces to a conventional Cattaneo- 
Mindlin solution with / = f s , so we have arbitrarily used this result to continue 
the curves to c = a [shown dotted]. 
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1.4 



Fig. 5: Comparison of the double Cattaneo-Mindlin solution with the JKR 
limit for i/) = 0.1. 

As predicted, the curve for A = 0.025 is very close to the JKR curve, though 
the maximum Q x is shifted slightly to the left. Notice incidentally that we 
might have chosen to plot the double Cattaneo-Mindlin curves as functions 
of the location of the mid-point (c + b)/2a of the slip-stick transition region, 
in which case this shift would be much reduced. For larger A, the maximum 
occurs at significantly lower values of c, but Qmax is still very well predicted by 
the JKR theory even for A = 0.4. 

Experimental measurements of static friction ceofficient are usually obtained 
by increasing the applied tangential force until sliding commences. However, 
it is clear that under these circumstances, microslip is likely to occur before 
gross sliding commences, and hence in the present geometry such experiments 
would lead to the static coefficient of friction being identified as Qm ax /P, which 
generally differs from f s . 
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1 


1.5 


2 


2.5 


3 



Fig. 6: The coefficient ratio f s /fd as a function of the apparent ratio Q max / fd P- 

Fig. 6 shows the relationship between f s /fd and the ‘apparent’ value of this 
ratio determined as Qmax/ fd P, for various values of 



(43) 


The dashed lines in this figure correspond to ranges in which the small-scale 
transition criterion s 0 c is n ot satisfied. We notice that the apparent static 
friction coefficient is always significantly lower than f s . The reason of course 
is that by the time Q max is reached, a significant part of the contact area has 
slipped sufficiently to transition to a local coefficient fd , and the measured 
coefficient is a weighted average over the whole contact area. 

Notice that the limiting case x = 0 can arise only if A = 0, meaning that the 
transition from f s to fd occurs over an infinitesimal slip distance. As explained 
in §3, the partial slip solution is then identical to the conventional Cattaneo- 
Mindlin solution with f = fd and hence slip occurs for Q = fdP regardless of 
the static coefficient of friction f s . This case is defined by the vertical axis in 
Fig. 6. 

7 Conclusions 

We have shown that the use of a slip-weakening friction law has a qualitative 
effect on the solution of microslip problems. The mechanics of the classical 
Cattaneo-Mindlin problem then have a mathematical structure similar to that 
of the adhesive contact problem, and we can identify an analogue of the ‘JKR’ 
solution, in which the extent of the stick zone is governed by the occurrence of a 
pressure-dependent mode II stress-intensity factor at the stick-slip boundary. By 
exploring the two-dimensional Hertzian geometry, we were able to identify the 
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equivalent fracture toughness, which is independent of the detailed goemetry, 
but proportional to the square root of the local contact pressure. We also defined 
a length scale So analogous to the small-scale yielding criterion whose value 
enables us to judge whether the singular solution gives a good approximation 
to the more exact solution. 

The tangential force reaches a maximum before the stick zone has shrunk to 
zero, at which point there will be a discontinuous change of state to gross sliding. 
This implies that estimates of the static coefficient of friction from experiments 
on the inception of sliding will generally significantly underestimate the values 
appropriate at the microscale. 
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